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ABSTRACT 

We examine the calculation of the solar neutrino flux based on the standard 
solar model (SSM). It is found that the solar neutrino data ( KAMIOKANDE 
experiment ) can be well described by the SSM with careful employment of nuclear 
data of Be{p,'~f)^B. The main point is that the simple-minded product ansatz 
of Coulomb plus nuclear parts should have a few percent uncertainties which 
induce the large reduction of the neutrino flux from ^B. Also, if the electron 
capture of ^Be inside the sun is suppressed, then the GALLEX experiment can 
be understood by the SSM calculation. 
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1. Introduction 

The solar neutrino problem is a long standing puzzle. The discrepancy 
between theoretical predictions of the neutrino flux by the standard solar 
model (SSM) and the observed data is still believed to be a factor of 2 
or more [1]. This problem, however, has produced many different kinds of 
refinements of solar internal structure model as well as new ideas in neutrino 
physics such as neutrino oscillations [2]. 

In this paper, we reexamine the calculation of the standard solar model by 
carefully considering the nuclear reaction data. To this claim, we may face 
criticisms that the nuclear reaction parts must have already been examined 
very carefully by all of the previous calculations. This is certainly right. 
The nuclear reaction data have been improved a lot and only those refined 
data have been employed. 

However, there is one important point which is required to reconsider in the 
previous calculations. That is, the Coulomb part calculated by the WKB 

method. The Coulomb coefficients can be calculated quite reliably if it is 
only one body problem. However, if it involves many body nature in the 
nuclear reaction, it is not very clear to what accuracy one can believe the 
WKB results even though we know that they cannot be very bad. 

Also, one knows in nuclear physics that the Coulomb problem is not as 
simple as one at first thinks. The Nolen-Schiffer anomaly is a good example 
[3] . The Coulomb displacement energy is not well described if one wants to 
discuss it to a very high accuracy [4]. 

Here, the problem is that the high energy neutrino flux is very sensitive to 
the Coulomb coefficients. In fact, the few percent change of the Coulomb 
coefficients may sometimes induce a large effect on the neutrino flux, leaving 
most of the solar structure quantities unchanged. In particular, the nuclear 
reaction data of ^Se(p, 7)^^ is most sensitive to the high energy part of 
the solar neutrino flux. As we will see below, a few percent increase of the 
Coulomb coefficient is enough to reduce the neutrino flux by a factor of 5. 
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Furthermore, the choice of the new Coulomb coefficient is perfectly consis- 
tent with the existing reaction data of '^Be{p,^yB [5]. At the same time, 
we reproduce all of the physical quantities of the solar internal structure at 
the same level of accuracy as the previous calculations of SSM. 

To summarize our results of the neutrino flux, we obtain the following 
neutrino capture rates for GALLEX [6], KAMIOKANDE [7], SAGE [8] 
and Homestake (Davis et al. [9] ). Here, BP95 and DS96 denote the recent 
calculations by Bahcall and Pinsonneaul [10], and Dar and Shaviv [11], 
respectively. We present the two different calculations ( Case I and Case II 
) which will be explained below in detail. 



Neutrino Flux 





Present cal. 


BP95 


DS96 


Experiment 




I 


II 


Homestake(SNU) 


4.5 


3.4 


9.3 ± 1.4 


4.1 ±1.2 


2.55 ±0.17 ±0.18 


KAMIOKANDE 
(lO^cm-^sec-i) 


2.9 


1.9 


6.62 


2.49 


2.73±0.17±0.34 


GALLEX(SNU) 


116 


114 


137±8 


115±6 


77.1±8.5l^:| 


SAGE(SNU) 


116 


114 


137±8 


115±6 


69±10l| 



2. The Standard Solar Model 

The internal structure of the sun is by now described reasonably well by the 
standard solar model. The chain of nuclear reactions is well understood. 
The description of the sun reduces to several couples of differential equa- 
tions which should be solved mostly by numerical calculations. Among the 
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parameters that enter in the equations, the opacity coefficient must be most 
ambiguous. However, recent studies to refine the SSM enable us to remove 
the ambiguity of the opacity fairly weU. In connection with the solar neu- 
trino problems, the ambiguity of the opacity may lead to a correction of a 
few tens of percents to the neutrino flux. In this respect, we have only a 
very little freedom left for neutrino flux [1]. 

The energy of the sun is governed by the nuclear reaction cross sections. 
The energy production rate ei2 for 1 + 2 — > 3 + 4 + Q reaction is described 

by 

QN1N2 <(TV> . 

ei2 = , . X erg/g • s 

(1 + 612) P 

where Ni and N2 are the number of particles in the reactions. 

Further, the < av > can be parametrized for nuclear reactions in the fol- 
lowing way except ^Be{e, VeYLi reaction. 



Na < av >= CiTg ^ exp f-CaTg ^ - f 



X (1 + C3T/ + C4TI + C5T9 + CsTI + CtTI) + CgTg ^ exp {-C.T,' 

(2.1) 

For the '^Be{e, v^YLi reaction, we employ the following form. 



Na < av >= 1.34 x 10^°To ^ 



X (^1 - 0.537r| + 3.867^ + I.2T9 + 0.0027r9 exp ^ O-O^^IS ^^ 2) 

Here, Na denotes Avogadro number. Tg is measured by 10^ K. 

The values of the parameters Ci, C9, Tq are determined from the nuclear 
reaction data and are listed in ref.[l,12]. 

The temperature T and the density p of the sun are determined by solving 
the following coupled equations. 
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dP _ GMp 
dr 

dM . 2 



dr 
dL 



dT 3KpL 



dr Wnacr^T^ 
dT 1 TdP 



(2.3a) 



= 47rr> (2.3&) 



^ = Anr^pe (2.3c) 
dr 



for radiative (2.3ci) 



for convective (2.3e) 



dr {n + l)ad P dr 

where P, M and L denote the pressure, the interior mass and the luminosity, 

respectively. Also, a and c arc radiation density constant and the velocity of 
light. Further, n and (n+ l)ad denote the opacity and adiabatic coefficient, 
respectively. 



3. ^Be{p, j^B reaction 

Now, we want to discuss the nuclear reaction '^Be{p, 7)^-6 since this is ob- 
viously the most important reaction to produce high energy neutrinos. In 
particular, we want to focus on the penetration factor Pcoui- This is ex- 
pressed in terms of the WKB calculation as 

C 

Pcoui = exp (-27r77) (3.1) 
where Co is a constant and rj can be described as 
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where denotes the reduced mass of the interacting particles. To the 
first order approximation, one may assume that the cross section can be 
described as a product of Pcoui and the nuclear part PnucI that is connected 
to the probability to make nuclear reactions. 

a{E) = -PcouiPnucI = exp{-27rr]) (3.3) 

where S{E) is a nuclear spectroscopic factor. 

Now, the question is to what accuracy we can believe the product ansatz 
of cq.(3.3) even though the WKB estimation is taken to be reliable. This 
is connected to the fact that the nuclear reaction of ''^Be{p,'~f)^B should be 
treated as a many body problem. Recent calculations by Brown et al.[13] 
show rather a large value of <S'(0). 

On the other hand, Xu et al. [14] claim that the S{0) value extracted 
from — > p + Be decay vertex constant is consistent with the observed 
value of Filippone et al. {S{Q) ~ 17.5 eV ) [5]. Thus, it is still far beyond 
determining the Coulomb coefficient to the accuracy of a few percent in the 
realistic nuclear many body calculations. 

Here, we do not want to rely on the simple-minded product ansatz of 
eq.(3.3). Instead, we assume the following form for a{E), 

(3.4) 

where A and B are free parameters which should be determined by repro- 
ducing the nuclear reaction data. We stress that our aim is not to reproduce 
theoretically the cross section data, but to find out some parameter sets that 
reproduce the observed cross section [5] . 

In fig.l, we show the comparison of the observed cross section of ^Be{p, ^)^B 
with that reproduced by eq.(3.4) with two choices of the parameter set A 



a{E) 



B 
E 



exp 



A 



^/E 
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and B. The first case (Case I) is the best fit to the nuclear cross section 
with the fixed value of A which is estimated by the WKB method. They 
are A = 4.70 x 10~^ erg^ and B = 17.5 eV. In the second case (Case II), 
we make the best fit to the nuclear cross section varying the values of the 
parameters A and B freely. We find that the best fit values of A and B are 
A = 4.80 X 10-2 erg5 and B = 20.0 eV. 

As can be seen, there are obviously some ambiguities which arises from 
the difficulty of the Coulomb cross sections once we want to understand it 
to a very high accuracy. With these two cases of the parameters, we can 
calculate the neutrino flux in the sun. 



4. The solar structure 

In the previous section, we have determined the parameters of the cross 
section < av > for the '^Be{p,'y)^B reaction. For other reaction cross 
sections, we have used the same values of parameters as those used in the 
calculation of Bahcall et al [15]. 

Here, we want to show our calculated result of the solar structure quantities. 
In fig. 2, we show the luminosity and the temperature of the sun as the 
function of the solar radius. The solid fines are the calculated results where 
the reaction cross section of ^i?e(j9, '~f )^B is used with the parameters A and 
B ( Case I ) as determined above. All the other nuclear data are the same 
as those used in the calculations of Bahcall et al. On the other hand, the 
dashed lines indicate the calculated luminosity and temperature by Bahcall 
et al. As can be seen from these figures, the shape of the luminosity and 
the temperature are almost the same between the two calculations. 

Therefore, we can conclude that the solar structure quantities are not so 
much influenced by the change of nuclear reaction data of Be{p,'y)^B, as 
expected. 
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5. The neutrino flux 

Since we know now how many reactions occur inside the sun, we can cal- 
culate the neutrino flux. 

In table 1, we show the neutrino fluxes as well as the capture rates at the 
Earth for GALLEX, KAMIOKANDE and Homestake experiments. In table 
la, we show the calculated results by Bahcall et al. [15] while, in table lb, 
the calculations by Dar and Shaviv [11] are shown. In table Ic, we show 
our calculated results with the Case I while, in table Id, the results with 
the Case II are shown. 

As can be seen from the table 1, the present calculations with the Case I are 
very similar to the ones by Dar and Shaviv. Therefore, it is confirmed that 
the KAMIOKANDE experiment is indeed consistent with the SSM calcu- 
lations with the careful employment of the nuclear reaction cross section of 

Further, the case II indicates that the ambiguity of the coulomb coefficient 
is so large that one has to be very careful for drawing any conclusions on 
the solar neutrino problems. At least, the result of the case II suggests that, 
once the ^Be neutrino flux is suppressed, then there is a fairly good chance 
that all the neutrino experiments fall into the range of the SSM predictions. 

For the Case II, one sees that the cross section of Be{p^ l)^B is best fitted. 
Here, the Coulomb coefficient is slightly different from the WKB value. In 
this parameter set, we find that the neutrino flux for KAMIOKAMDE is 
a little bit too small compared to the data. Instead, the Homestake and 
GALLEX experiments will be in the range of the present calculation once 
the ^ Be neutrino flux is suppressed. 

Also, in the Case II, the 5(0) value is found to be 5(0) = 20 eV . This 
suggests that the 5(0) factor depends on the factorization ansatz of eq.(3.3). 
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6. Conclusions 

We have presented a new calculation of the standard solar model with the 
emphasis on the careful considerations of the nuclear reactions of '^Be{p, ■j)^B 

We show here that the solar neutrino capture rates are consistent with 
the observed data for the KAMIOKANDE experiments. We believe that 
possible refinements may improve the accuracy of the neutrino capture rates 
by 20 ~ 30 % so that the GALLEX experiments may well be in the range 
of the SSM picture. In particular, the suppression of the ^Be electron 
capture inside the sun will lead to the understanding of the GALLEX and 
Homestake experiments in a natural way. 

Therefore, we conclude that the solar neutrino fluxes are mostly consistent 
with the standard solar model with careful considerations of the nuclear 
reactions of ^Be{p, 7)^-6. 

In the course of the present study, we received a preprint of the new cal- 
culation by Dar and Shaviv which shows very similar results to the present 
calculations. This confirms that the present result does not so much de- 
pend on the modeling of the sun as far as we take into account the gross 
structure of the sun. 
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Table captions 

Table 1 : We plot the calculated neutrino flux from various nuclear reactions 
together with the experiments. Table la shows the calculation by Bahcall 
and Ulrich [15] while Table lb plots the calculation by Dar and Shaviv [11]. 
Tables Ic and Id are the present calculations with the parameter sets of 
Case I and Case II, respectively. 



Figure captions 

Fig.l : We show the nuclear cross section of ^Be{p, 7)^-6. The black circles with 
error bars are the observed data by Fihppone et al [5] . The solid line is our 
calculation with the Case I parameters while the dashed hne with the Case 
II parameters. 

Fig. 2 : The properties of the internal structure of the sun are shown as the 
function of the radius. The solid lines show the present calculations while 
the dashed hues the ones by Bahcall et al [15]. The L, M, T, p and P denote 
the luminosity, the mass, the temperature, the density and the pressure of 
the sun. 
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Table la 



source 


Flux 

(cm~^sec~M 


Homestake 


GALLEX 

1 kJi >l KJ J 


BP88 

SAGE KAMIOKANDE 
(SNU) (lO^cm-^sec-i) 


pp 


6.0 X 10^° 




70.8 


70.8 




pep 


1.4 X 10^ 


0.2 


3.0 


3.0 




^Be 


4.7 X 10^ 


1.1 


34.3 


34.3 




8B 


5.8 X 10^ 


6.1 


14.0 


14.0 


5.8 


13N 


6.1 X 10^ 


0.1 


3.8 


3.8 




150 


5.2 X 10^ 


0.3 


6.1 


6.1 




Total 




7.9 


132 


132 


5.8 


Experiment 




2.55 ±0.25 


77.1 ± 8.5l|:^ 


69±10lf 2.73 ±0.17 ±0.34 



Table lb 



source 


Flux 
(cm~^sec~^) 


Homestake 
(SNU) 


GALLEX 
(SNU) 


DS96 

SAGE KAMIOKANDE 
(SNU) (lO^cm-^sec-i) 


VP 


6.1 X 10^° 




72.0 


72.0 




pep 


1.43 X 10^ 


0.20 


3.06 


3.06 




^Be 


3.71 X 10^ 


0.87 


27.1 


27.1 




«B 


2.49 X 10^ 


2.62 


6.01 


6.01 


2.49 


13N 


3.82 X 10^ 


0.06 


2.38 


2.38 




150 


3.74 X 10« 


0.22 


4.39 


4.39 




Total 




4.1 


115 


115 


2.49 


Experiment 




2.55 ±0.25 


77.1±8.5l|^ 


69±10lf 2.73 ±0.17 ±0.34 
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Table Ic 



source 


Flux 

(cm~^sec~M 


Homestake 


GALLEX SAGE KAMIOKANDE 
(SNU) (SNU) (lO^cm-^sec-i) 


pp 


5.7 X 10^° 




67.3 


67.3 




pep 


1.4 X 10^ 


0.20 


3.00 


3.00 




^Be 


4.7 X 10^ 


1.10 


34.3 


34.3 




8B 


2.9 X 10^ 


3.05 


7.00 


7.00 


2.9 


13N 


3.7 X 10^ 


0.06 


2.30 


2.30 




150 


2.2 X 10^ 


0.13 


2.58 


2.58 




Total 




4.5 


116 


116 


2.9 


Experiment 




2.55 ±0.25 


77.1 ± 8.5l|:^ 


69±10lf 2.73 ±0.17 ±0.34 



Table Id 



source 


Flux 
(cm~^sec~^) 


Homestake 
(SNU) 


GALLEX SAGE KAMIOKANDE 
(SNU) (SNU) (lO^cm-^sec-i) 


VP 


5.7 X 10^° 




67.3 


67.3 




pep 


1.4 X 10^ 


0.20 


3.00 


3.00 




^Be 


4.7 X 10^ 


1.10 


34.3 


34.3 




«B 


1.9 X 10^ 


1.95 


4.47 


4.47 


1.9 


13N 


3.7 X 10« 


0.06 


2.30 


2.30 




150 


2.2 X 10^ 


0.13 


2.58 


2.58 




Total 




3.4 


114 


114 


1.9 


Experiment 




2.55 ±0.25 


77.1±8.5l|^ 


69±10lf 2.73 ±0.17 ±0.34 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9707092vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9707092vl 



